Note on MacLennan-Zubarev Ensembles and QuasiStatic Processes 



Shuichi Tasaki a 

Taku Matsui 6 

"■Department of Applied Physics, Waseda University, 
Graduate School of Mathematics, Kyushu University 

Relation between natural nonequilibrium steady states of Ruelle and MacLennan-Zubarev 
. ensembles is discussed for a system consisting of M infinitely extended systems coupled with 

a small system. And a characterization of quasistatic processes is investigated for a small 
system coupled with a single reservoir. 

>v 



in 



43 



INTRODUCTION 



The understanding of irreversible phenomena including nonequilibrium steady states is a long- 
standing problem of statistical mechanics. Recent progress in the research of mesoscopic systems 
brings a new aspect into this problem. In these systems, coherence (a quantum dynamical aspect) 
may be observed in a dissipative transport (an irreversible phenomenon) and the two aspects should 
be discussed simultaneously. Usually, a mesoscopic system couples with much larger environments 
and the interaction is not weak. As a result, the system cannot be clearly distinguished from the 
environments. Therefore, it is natural to deal with a mesoscopic system plus its environments as 
an infinitely extended system. 

Statistical mechanics of infinitely extended systems has been developed so far 0, Q, 0| and, 
recently, their nonequilibrium properties are studied intensively. Those include analytical studies 
on nonequilibrium steady states of harmonic crystals [J, E| , a one-dimensional gas [f| , unharmonic 
chains I, an isotropic XY-chain g, systems with asymptotic abelianness g , a one-dimensional 
quantum conductor [icj|. an interacting fermion-spin system [ill ], fermionic junction systems [l2T ]. 
a quasi-spin model of superconductors and a bosonic junction system with Bose-Einstein 
condensates 0]. Entropy production has been rigorously studied as well (see Refs. P. Hll. Hil ull 



171 la ] , and the references therein) . See also reviews [19|, |2fJ . 
^ ■ Fta C*-dynamical systems with ^-asymptotic abdlan properties, Ruellei showed the exis- 

tence and stability of nonequlibrium steady states (NESS) which are naturally obtained from lo- 
cal equilibrium ensembles. We previously announced |2ol| that, under stronger conditions, these 
natural NESS can be regarded as nonequilibrium ensembles proposed by MacLennan[2l| and 
Zubaxev0. In this article, we give the proof of this statement. On the other hand, one of 
us (ST) investigated 23] a thermodynamic behavior of a driven small system coupled with an 
infinitely extended reservoir. Another purpose of this article is to give the detail proof and general- 
ization of this observation. The rest part is arranged as follows: In the next section, we summarize 
the setting. In Sec. IIII1 the relation between natural NESS and Maclennan-Zubarev ensembles is 
discussed. Sec. IIVI is devoted to the discussions of the thermodynamic bahavior of a small system 
coupled with an infinitely extended reservoir. Some remarks are given in Sec. 
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II. C*-DYNAMICAL SYSTEMS 



A. Field algebra 



The system <S in question is described by a collection J- of all finite observables called a C*- 
algebra, which is a complete linear space with a norm || • ||, where a procduct AB and antilinear 
involution *: A — * A* (\/A,B G J~) are defined and whose norm satisfies < ||^4||||5|| and 

the C*-property: = ||^4|| 2 . We consider the case where an autonomous time evolution is 

defined which is described by a strongly continuous one-parameter group of *-automorphisms Tt 
(t G R). Namely, Tt is a linear map satisfying T t (AB) = T t (A)Tt(B), T t (A*) = r t (A)* , tq = I (I: the 
identity map), T t r s = T t + S and lim^o ||r((^4) — A\\ = (^A G J 7 ). Then, according to the theory of 
semigroups, there exists a densely defined generator 5 of 



lim 



6(A) - l{ n (A) - A} 



, (Vie D(6)) 



where D(5) is the domain of 5. 

We further assume that two more *-automorphisms are defined: 



(i) a strongly continuous L-parameter group of *-automorphisms ((p G R L ) satisfying 
tp 1 + l fi 2 ' wri i cri represents the gauge transformation. 



(ii) an involutive ""-automorphism which is represented as & = a^ with some tpo G R 



L 



The automorphisms Tt, and G are assumed to commute with each other: 
@r t = r t 6, Qct0 = a^Q, T t a0 = a^r t ( y t G R , V <p G R L ) . 

A subalgebra A C J- consisting of invariant elements under the gauge transformations 
(Cp G R L ) is called the observable algebra, which describes observable physical quantities. The 
*-automorphism defines the even and odd subalgebras, respectively, F + and JF_: 

J-± = {A G T; e(A) = ±A} . 

When the system involves fermions, even and odd subalgebras correspond to dynamical variables 
which are sums of products of, respectively, even and odd numbers of fermion creation and/or 
annihilation operators. Note that, since a s g x (s G R, ey. a unit vector whose Ath element is unity) 
defines a strongly continuous group of *-automorphisms, it has a densely defined generator which 
is denoted as g\. The C*-algebra T with these *-automorphisms is called a field algebra 0, H3j- 



B. Decomposition of the system and initial local equilibrium states 

We consider the situation where the system <S can be decomposed into M independent infinitely 
extended subsystems TZj (j = 1, ■••M), which play a role of reservoirs, and a finite-degree-of- 
freedom subsystem So interacting with all the others. More precisely, the algebra J 7 is a tensor 
product of M infinite dimensional subalgebras Tj (j = 1, • • • M), corresponding to TZj, and a finite 
dimensional subalgebra J-s, corresponding to Sq: 

F = F s ®J : i®---®J : M , (1) 
such that the following conditions are satisfied: 
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(SI) There exists a gauge-invariant time evolution group r t (t G R) which is a perturbation of 
r t by a selfadjoint element — V G A PI -D(<5) and which is a product of strongly continuous 

groups (J = 1, • • ■ M) independently acting on subalgebras J-j (j = 1, • • • M) 

,(1) _(Af) 



(2) 

t:[S} 



(7) 

Namely, f} s commute with each other and leave the elements of other subalgebras 
invariant:^ 

f t (i) (A) = A ( V A G fc ^ j, k,j = l,--- M) (3) 

f .U) f {k) = ~(k)~(j) (t,s £~R, k ^ j, k,j = 1, ■ ■ ■ M) . (4) 

(52) The gauge *-automorphism is a product of strongly continuous groups 6$ 
(j = !,-■■ M) and d^: 

= <x$aty • • • Of > ( 5 ) 
where a9 and a? independently act on subalgebras Tj (j = 1, • • • M) and Ts 

af(A) = A ( V A € J^, fc ^ j, k,j = S,l,--- M) (6) 

= ag>ag (<^ 2 GR L , Mi, fe,i = 5, 1, • ■ ■ M) (7) 

(7) (fc) 

And the groups f f and dV commute with each other: 

f^af = afr^ (j, k = 1, • • • , M; i G R, £ G R L ) . 

States are introduced by listing expectation values. Namely, each state is identified with a 
linear map oj from i G J to an expectation value u){A). The positivity condition uj{A*A) > 
and normalization condition w(l) = 1 (with 1 G .T 7 the unity) are required. As in the previous 
works |4l. la. lid llll. Ila. lla. Il7| . we are interested in the evolution of initial states where M infinitely 
extended subsystems (reservoirs) are in equilibrium with different temperatures and different chem- 
ical potentials and the finite subsystem is in an arbitrary non-singular state. As discussed in [Til . 17], 
such states are specified as a KMS state: 

(53) Let a% {x € R) be a strongly continuous group defined by 

M 

i=i 

where (5j and flj = ((Jj , ■ ■ ■ are, respectively, the inverse temperature and the set of 
chemical potentials of the jth reservoir. The operator D$ (G J-s H A) is selfadjoint. Then 
an initial state uj, which we are interested in, is a KMS state with temperature —1 with 
respect to a^. Namely, uj is a state such that, for any pair A,BeJ-, there exists a function 
Fa,b{x) of x-analytic in the stripe {x G C; > Imx > —1} and satisfies the KMS boundary 
condition: 

F A)B {x) = u(Aa%{B)) F A! b(x - i) = uj(<t%{B)A) (x G R) . (9) 
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(i) (i) 

Hereafter, we assume that the domains of the generators Sj and gV , respectively, of f t and 
a^g (t, s € R) contain the domain D(5) of the generator 5 of r t : 



(S4) D(5) C D(5j) D(5) C D<$) ( v j = 0, 1, • • • M , V A = 1, • • • L). 

Then, the condition (SI) implies that the domain of the generator 5 V of tY is identical with D(5): 
D(5 V ) = D(5) and 

5(A) = 8 v {A)+i[V,A] (A e D(6)) (10) 

N 

5 v {A)=Y.h(*) {A e D(S)) (11) 

Moreover, the domain of the generator 5^ of o% also includes D(5) and 

N 

6 W (A) = (s j (A)-^ ) ~g?(A))}+ i [D s ,A} . (A 6 D(S)) (12) 

i=i 

Note that a decomposition without the finite subsystem is possible as well. 



C. i 1 -asymptotic abelian property 

In thermodynamics, environments are assumed to stay in equilibrium under arbitrary processes 
and their details are considered to be unimportant. Hence, thermodynamic environments would 
be well-modelled by systems with appropriate ergodicity. As one of such an example, we consider 
systems satisfying the L 1 -asymptotic abelian property. 

The time evolution is said to satisfy the L 1 (C?)-asymptotic abelian property if there exists a 
norm dense *-subalgebra Q such that 

»+oo 



/+oo r+oo 
dt\\[A,T t {B)}\\= / dt\\Ar t {B) - T t {B)A\\ < +oo (A G Q, B G Q n F+) 

/+oo r+oo 
dt\\[A,T t (B)]+\\= dt\\AT t {B) + T t {B)A\\ < +oo (A, B e G n 
-oo J —oo 

where T± are even/odd subalgebras. This property implies rapid decay of correlations and is 
satisfied by free fermions in R d (d > 1) (Example 5.4.9 of Ref. 0]). Note that, if a system 
admits bound states, it does not satisfy the L 1 -asymptotic abelian condition as there exist bounded 
constants of motion, i.e., observables C satisfying T t (C) = C. 

III. NESS AND MACLENNAN-ZUBAREV ENSEMBLES 

The state Ut at time t starting from the initial state u is given by 

u t (A)=u(r t (A)) , ?AeF) (13) 

and, under the setting (S1)-(S4), nonequlibrium steady states are expected to be obtained as its 
weak limits for t — > ±oo. As shown by Ruelle[9j, it is indeed the case. Namely, under the setting 



5 



(S1)-(S4), if the time evolution n is L 1 (<7)-asymptotic abelian and the perturbation V is an element 
of Q, the limits 

lim u (T t {A)) = w(7±(A)) = u±(A) 

t— >±oo 

exist for all A € and define nonequilibrium steady states, where j± are M0ller morphisms defined 
by 



lim 

t— >±oo 



V t (A)- 7± (A) 







Also Ruelle showed the independence of the limits on the way of separating the whole system into 
reservoirs and a small systemj^]. See also Refs. [ill. In this section, as announced in Ref. |2ol |. 
under a stronger condition, these steady states ui± are shown to be inerpreted as MacLennan- 
Zubarev ensembles 21, zl. Indeed, one has: 



Proposition 1: KMS characterization of evolving states 

Under the setting (S1)-(S4), the state ojt at time t is a KMS state at temperature —1 
with respect to the strongly continuous group of *-automorphisms 



(14) 



where jt = Tj, and its generator is given by 



S®(A) = 6 u (A)+i I ds 



A 



( v AeD(^) = D(5 UJ )) 



(15) 



Proposition 2: KMS characterization of steady states 

Under the setting (S1)-(S4), if the time evolution *-automorphism r% is L l (Q)- 
asymptotically abelian, V £ Q and the M0ller morphisims j± are invertible, the steady 
states u± are KMS states at temperature —1 with respect to the strongly continuous 
group of *-automorphisms 



(16) 



Furthermore, if d~uj(V) € Q, its generator 5^ satisfies 



S±(A)=6 u (A) + i 



ds 



=Foo 



(Uv) 



A 



1 A g d(5 w ) n g) 



(17) 



NB 3 For finite systems, the KMS state u with respect to the *-automorphism corresponds 
to the density matrix 

1 N L 

^4°p{-E^-£W A) )}' 

3=1 A=l 

where Z is the normalization constant, [3j, Hj, fj,^ and are, respectively, the local tempera- 
ture, local energy, local chemical potential and local number operator of the jth reservoir. Then, 
as a result of the Liouville-von Neumann eqation, the density matrix at time t is given by 

-. AI L 

3=1 X=l 
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-.ML .q 

1 exp{- £ &i [Hi ~ E M?W A) " / ^ S (J/)] } 
i=i a=i 17 ~* 



where Jj = |r t (^) - E^i ^^(Jvf ) 



-i , V] + Ea=i Wj X) , V] is non 



t=o 



systematic energy flow, or heat flow, to the jth reservoir and we have used 



JV. 



(A) 



t^Hj) = Hj - ds-r s {H 3 ) , r_ f (Af } ) . , 

On the other hand, as <5j(V") = i[Hj, V] and = i[Nj X \ V], one has 
M M L 

E/V^D = -r.(5>{*i(V) - £/W(v)}) = -r. (* w (y) 

i=i i=i 



and, thus, 



A=l 



-, AT L 

= \ exp{- £ ft - £ /^ivf + 



A=l 



which is a KMS state generated by ^(A) + i J ds[r s (Su(V)^ , A] . This observation is nothing 
but the finite dimensional version of Proposition 1. 

NB 4 For infinite systems, an interesting case is that where the right-hand side of (|17|) generates 
cr^. Then, if the integral 



=poo 



dsr s (6UV) 



(18) 



would converge, uj± would be perturbed KMS states of the initial state u by self-adjoint operators 
V±. Moreover, NB 3 suggests that the corresponding density matrices would be 



P± = 7 «xp{- E Pi ~ E M?^ A) - / dsr s (J])] } . 

j=l A=l 17=1:00 



(19) 



Note that such ensembles for steady states were introduced by MacLennan[2l|| and Zubarev[2^ 

However, if the steady states carry nonvanishing entropy production, the integral V± does not 
converge since the steady-state average of its integrand is nothing but the nonvanishing enropy 
production rate at the steady states: 



M 



fim oj(t s (S w (yj)) = w±(E 13, J]) / . 

i=l 



This observation is consistent with the results by Jaksic and Pilletjl7l| who showed that, if the 
steady-state entropy production is nonvanishing, the steady state is 'singular' with respect to the 
initial local equilibrium state. Thus, the original proposal (|19|) by MacLennan and Zubarev cannot 
be justified. Rather, the KMS states with respect to generated by (|T7|) should be regarded as 
a precise definition of the MacLennan-Zubarev ensembles. 

Proof of Proposition 1 
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As easily seen from (SI) and (S4), the initial state u is invariant under r t . And one has 
u>t(A) = lo (jt(A)) = lo ('yt(A)) ( V A G T) ■ Since 74 is a *-automorphism and u is a KMS state 
with respect to at temperature —1, ut is a KMS state with respect to ■y^ 1 a l ^-yt at temperature 
-1 (cf. Prop. 5.3.33 of Ref . 0). 

Now we consider the generator. In terms of the one-parameter family Yj of unitary elements 
defined as a norm convergent series: 



n=l 



dt n T. tn (V)---T_ t2 (V)T_ tl (V) 



(20) 



one has jt(A) = Y t AY t *. Thus, jt(A) A G T) is differentiable with respect to t and —^t{A) 
*7t ([T_ t (V), A]) . This leads to 

d{ 7t -V£ 7 ^)} 



ibr^Ttr-tCV) - r_ t (y), 74-^7*^)] = * ^ - V) ,7^7^) 



dt 



or 



7r i ^7t(^) = + i / - V) nl^lsiA) 



(21) 



where 7fT_ t = rY t , T V t a% = &x T -t anc ^ 7t lr -t = T -t have been used. Since V G D(<5) C D(5u), 
one has the following relation for any A G D{5 UJ ): 



lim 



x=0 



.r 



+ if ds 



7s ^x7«(^J 



$ w (A) + i / ds [T a [5 w (V)),A 



(22) 



This suggests that 5 wt defined by 



S m (A) = L(A)+iJ ds [t s (L(V)) , A 

is the generator 6$ of ^^o^^t- As we see, it is the case. First, we note that, since f® t dsr s {5 U {V 

is selfadjoint, 5 Wt generates a strongly continuous group (cf. Theorem 4.1 of Ref. 3] ). Then as a 
result of the theory of semigroups (Prop. 3.1.6 of Ref. ji| ), its resolvent {fil + ^ t ) _1 (/•* G R-\{0}) 
is bounded, where 1 stands for the identity operator on T. For any B G J 7 , let A = (fil + 5^ t )~ l B, 
then i e and (J22J) leads to 



( M l + ^)A = M+^7r 1 ^7*(^) 



5 



x=0 

-1 



On the other hand, as a generator of strongly continuous group j t L <?xlt, the resolvent (fil + sffl) 1 
of 6$ is again bounded. Therefore, one finally has 

{pi + V$y x B = A=(jil + 5 Ut y l B 
for all B G ^ and <5$ = 4, t . Moreover, as £>(<L t ) = = D(6 U ). 
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Proof of Proposition 2 

Remind that ui±(A) = uj (^±(A)) y 'A G T). When the M0ller morphisms 7± are invertible, 
they are *-automorphisms. Thus, since u> is a KMS state with respect to at temperature — 1, 
u± is a KMS state with respect to 7± 1 0"£7± at temperature —1 (cf. Prop. 5.3.33 of Ref. 

Now we show the second half. Since *-automorphisms are norm preserving (cf. e.g. Corollary 
2.3.4 of Ref. H ), one has 



7iX 7 ±n(A) -7± < 7 ±(A)|| < h± n (A)- 1± (A)\\+F(x,n) 



(23) 



where 7± n = T± n V± n and 



n) 

As the M0ller morphism is a strong operator limit of 7± n , 

lim [| 7 ±n(A)-7±(^)|| = 0, 

n— >+oo 

lim F(x, n) = . 

Next, result of 



7±n 7iX7±(^) "7± h±o%>y±(A) 



CA G ^) 
(for each x) 



(24) 
(25) 



|F(x,n)-F(x',n)| < || 7 ±„7i 1 K7±W - <#7±(^)} II 
= K 7 ±(A)-<7±(>1)|| , 

the function a; — > F(x, n) is continuous uniformly with respect to n. Hence, the limit (|25|) is 
uniform with respect to x on any finite interval in R. In short, 



lim || 7 ^< 7±n (A) - 7± ± < 7± (A)|| = 



(26) 



where the convergence is uniform with respect to x on any finite interval in R. Hence, as a result 
of Theorem 3.1.28 of Ref. [j] , the generator 5^ of 'y^ 1 a^j± is the graph limit of the generators 



of -y ±n a^j± n . Namely, if a sequence {A n }+^ with A n G D(5l ±n) ) satisfies A n -> A and 



(A 



B (n — > +00) , then one has -B = <5^(j4). 



For each A G D(5 U ) fl £7, let A n = A, then obviously linin^+oo = A, A n G -D(<5i ±n ^) because 
of Proposition 1. And, as ^(V) G Q is assumed, the L 1 (£/)-asymptotic abelian property of r t 
implies 



lim 5^ n \A n ) = 5 u) {A) + i 



=Foo 



(LOO), a 



Therefore, one has the desired result 



S±(A) = 6 u (A)+i 



ds 



=poo 



TslSUV)) ,A 



IV. QUASISTATIC PROCESS AND CLAUSIUS EQUALITY 

In this section, we consider the case where Ts is the algebra of bounded operators on a finite 
dimensional Hilbert space and the small system couples with a single reservoir, namely the case of 
M = 1. Then, the total system is described by the tensor product !Fs ® T\ and the 'decoupled' 
evolution tY acts only on the reservoir algebra: tY = T t . Since the system is finite dimensional, 
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(S) 

the generator of the gauge transformation a s g is a commutator with a self-adjoint element N\ (8> li 
where N\ G .Fs is the number operator of the Ath particles and li is the unit of T\. The system- 
reservoir interaction is assumed to be proportional to a coupling constant k: kV where V E D(5) 
is self-adjoint and gauge- invariant: a^iV) = V. Note that the gauge invariance of V leads to 

$\V) + i[N x ® li, V] = (A = 1, • • • L). 

We are interested in the response of the whole system under a time-dependent perturba- 
tion: W(t) (g) li where W(t) € !Fs is twice continuously differentiable in norm, W(t) = Wo for 
t < and [N\ <S> l\,W(t) ® li] = 0. Initially, the whole system is prepared to be an equilibrium 
state ijj of the inverse temperature j3 and the chemical potentials p\, namely, a KMS state at (3 

(i) (i) 

with respect to a s = a s a_ s ^ where 04 is defined by 



^dt^ = ^ + i[W ° + KV,A] ) ' a t\ A )\t=o=A, ( v AeD(5i)). 



The time evolution r t is given by the solution of 



(A) w (1 



dt 



TFfcW+iiKV + WWsiuA]), T^(A)\ t=0 = A, (A e D(5 X )) , 



and the state at time i by u^(^4) = w (t^ (A)^j A € JF). Now, we define the the system-energy 
increase Zt induced by the mass flow, the work Wt done on the system and the heat Qt absorbed 
by the system during the time interval T as follows: 

L 

Z T = X)ma{^(JVa®1i)-w(JVa®1i)} 

A=l 

W T = J dt cot (iW(t) ® U] 

Qt = {u>t(W(T) ® li + kV) - uj(W ® li + kF)} - W t - Z t 

Then, one has the following results: 

Proposition 5: Stepwise perturbation 

Suppose W(t) = Wf (t > to > 0). Let be the evolution defined by 

(/) 

^d/^ = "* /} + + kV ' A] ) ' a F ) (^)l*=o = ^, (Ae4)), 

wj- be a KMS state with respect to a s = ct\ ,a- 8 p at (3, and (TCf,iTf,Qf) be its 



GNS representation [3 ll| , Then, if the initial state u is the unique KMS state and a 

self-adjoint operator L f on Hf defined by e^/V/ 0*1)0/ = tt / (q;^ ) (vl € .F), 

which will be called the Liouvillian of a\ , has a simple eigenvalue at zero and the 
absolute continuous spectrum, one has 

lim lim 0Q T = S(p f ) - S(jh) - S(p f \p U) ) 

K — >U 1 — > + CsO 

where S(p v ) = — Trjp^ lnp^} [y = i, f) is the von Neumann entropy of the density 
matrix p u , S(pf\pt ) = Tr{p to (ln p to — Inpj)} > is the relative entropy between 
density matrices p f and p to , and pi = exp(-/3(W - Ea MA^Yx))/Si> Pt = u t Q Pi u t , 
Pf = exp(— /3(Wf — J2\ Ma^Ya))/S/, with Sj, 5/ the grand partition functions. The 
unitary element u t E Fs is the solution of ^u t = iu t W(t), u t \t=o = Is- 
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Proposition 6: Staircase perturbation 



Suppose that T = J2j=i Tj and the interaction W{t) has a staircase form: W(t) = 
W + (j-l+ip(t-f j _ 1 ))(W f -W )/N for %_ x <t<f 3 where % = ££ =1 T k , (f„ = 0) 
and tp(t) is a twice continuously differentiable real-valued function with tp(0) = 0, 

tp(t) = 1 for t > t . Define Wj = Wo + j{Wf — Wq)/N and introduce the group 
by 

da US - a U) r 6l{A) + i[w . + KVf A] \ ; a P(A)\ t=0 =A, CAeDfa)). 



dt 

Let ujj be a KMS state with respect to cr£^ = ai^a_ s ^ at /?, and (Hj,iTj,Qj) be 
its GNS representation. Then, if the initial state lu is the unique KMS state and 
the Liouvillian Lj defined by e iL ^irj(A)Slj = 7Tj(e4 (A))Slj (A G JF) has a simple 
eigenvalue at zero and the absolute continuous spectrum, one has 

2%A.~ = S(P,)-S(ft)+o(l 



NB 7: Because of the return-to-equilibrium property, the final state (and every intermediate state 
in Proposition 6) of the whole system is an equlibrium state. The limit of n — > implies that the 
coupling between the system and the reservoir is negligibly small. Thus, the processes treated in 
Proposition 5 and Proposition 6 precisely correspond to those in the classical thermodynamics |2(|. 
The Clausius inequality follows from Proposition 5: lim K ^o hm:r->+oo @Qt < S{Pf) ~ S{Pi) because 
of the positivity of the relative entropy. On the other hand, Proposition 6 implies that, if the whole 
system changes very slowly (N 3> 1) so that the whole system is in equlibrium at every instant, the 
Clausius equality holds. Thus, the thermodynamic entropy is given by the von Neumann entropy, 
as expected, and the process described in Proposition 6 is nothing but a quasistatic process. Note 
that, as kV is responsible for the equlibration, the weak coupling limit k — > should be taken after 
the long term limits. 

NB 8: The two propositions deal with the entropy change of a subsystem in contrast to the 
previous works 0, Till IliL llfil . 12 , 0, E3] , which have discussed the entropy production of the 
whole system. In this respect, the present work shares a common interest with that by Maes and 
Tasakip3|. who derived a relation equivalent to the Clausius inequality from dynamics for a class 
of finite classical systems. The difference between the present setting and that of Ref. .28] lies 
in the fact that we deal with an infinite-degree-of-freedom reservoir in order to have equilibrium 
states as the final state of Proposition 5 and as the intermediate states of Proposition 6. 

NB 9: Frohlich, Merkli, Schwarz, and Ueltschi[27| have discussed the Clausius equality for 'adi- 
abatic' processes where the difference between the true and the reference states |32j restricted to 
a certain subsystem is small for all times. However, conditions on the dynamics which realize the 
'adiabatic' processes are not given. Instead, Proposition 6 gives a concrete example of quasi-static 
processes, which is similar to the one studied for the classical stochastic systems by Sekimoto (2fll |. 

Proof of Proposition 5: The heat is rewritten and, then, the assertions are shown. 
Expression of heat: Let T t be solutions of -^T t = iT t f^ l \w{t) (B> li + kV), T t \t=o = 1, where 
1 = l s ® li, then rf(A) = Y t fp(A)Y* t and 

.d 



dt 



* - E rf)(r t )rr) = (5i - £ rfxr^Wt) ® i x + K v))T* t 

A=l A=l 



" E ^xg ( xWt)fr(W(t) ® l! + KV)T* t = r^((5i - E rfX^)) (27) 

A=l A=l 
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where we have used (Si - £a=i fix~g\ L) )(AB) = (<5i - Ea=1 »\9x)( a ) b + A (h ~ £a=i 
*i(W(t) ® li) = and ^ 1} (W(t) ® l x ) = 0. Reminding 

d 



d_ 

dt 



—r t w (N x ® li) = rf (<5i(iV A ® l x ) + i[W(t) ® li + kV, N x ® 1] 

= rf (t[«7,i^ ® li]) = rf (^(kF)) , (28) 

T t w (W(t) ®li + «V) = T^(s 1 (W(t)®l 1 +KV) + i[W(t)®li + KV,W(t)®l 1 + KV] 

1.) = rJ"(ft(«V)) 4- r-(«a 1.) , (29) 



one obtains 



d 

-i-r 
dt 



v A=l / 

drf(W(t) ® lj + kV) w (dW{t) \ ^ dl W( Nx % x 
A Tt l^^ 0ll J-^^ A A 

A — 1 



(30) 



and, thus, 



Q T = lo(t^(W(T) ® li + «V)) - w(W(0) ® li + «F) 



£ dsu; (rf ® li)) - f> A (w(rf (JV A ® li)) - u(N x ® l x ) 

A— 1 

'^((5i-Erf)(rr)r^ . (3i) 



Proo/ o/t/ie assertions: Let be the solution of ^r[ f) = «rpV t (1) (W/®li + rcF), r^ /} | t=0 = 1, 
then, as in the derivation of (|3U[). one obtains (A) = r^(.A)r£^*, 

{fr - E rf)^)^*} = 4 f) {(Si - Erfx^)) 

= ^a| /} ® li + kV- E>^ ® 1^ 



and, thus, 

L 



- - MAsfXrFVP = - E ^a^a) ® ii + kV) - (w f - E ^a^a) ® ii - «y 

A=l A=l A=l 

Then, because of Ft = Tt f^ (j^i-to) f° r * — ^o> the heat flow to the reservoir is rewritten as 

= w (r%f to a^l to ((w f - £ /^a^Va) ® ii + «^)r t * r«^ 
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— ioj 



{h-jz^ ] ){Tt )T* t \ 
v A=l / 



(32) 



where a unitary element 1^ is the solution of ^j+F^ 1 = iT\ l> f^ x> (Wq <8> li + kV) with T 



t |t=o — 1, 



I\ = fi 1 '* 1 (r^ ) is again unitary, and, in deriving the first two terms, we have used a relation 



u(ft (A)) = uj(a y t l> (f^(A) )) = ^(T^AT^'i) resulting from the a^-invariance of u> 



(0 (zMi 



Let us begin with the evaluation of the T-independent terms. Since 1(27)1 gives 

- - £ rfxr^jr* = K f° dsrY ((h - E mPm) , (33) 

A=l J0 V A=l / 



the last term is evaluated as 

L 



A=l / 



< 



A=l 



to 

< |/c| / ds 



k to 



(*i-Erf)(v) 



A=l 



(34) 



Since ||r io || = 1 and Ht^^V)!! = ||V||, the third term is bounded by |k|||V||: 



lo 



< \k\ \\V\\ . 



By taking into account the fact that Tf(W(t) <S> li) = W(t) ® li and comparing the equations 
for and Ut, one obtains 



r 4o = u to ® li + Ar to , Ar to = i« / ° dsr s r«(F) ) ® li) . 

J 



(35) 



Then, because of f^\(W f - Ea=i^A^a) ® li) = (W) - Ea=i^aNa) <8> li, the second term is 
evaluated as 



lo 



(T to f£\(W f - VxNx) ® li)r t * ) - u/(K,(W) - Ma^aK) ® li 



A=l 
L 



A=l 
L 



< 



LO 



LO 



T to((w f - J2 ^a^Va) ® ii)r* 

A=l 
L 

(ATtoCCWy-EMA^A)®!!)^ 



a; 



+ 



K(^/-E^a)<)®1i 



A=l 



A=l 
L 



lO 



«to(W>-EMAJVA))®liAr; 



A=l 
L 



< 2||Ar t0 || ||(W> - E Ma^a) ® liH < 2|K|t ||y|| - E Ma-^a) ® 1i|| , 

A=l A=l 

where we have used ||Ar t() || < |/c|io||V||. 

Now we turn to the first term, which is estimated with the aid of the following lemma: 
Lemma: 

lim u)(Aa[ f) (B)C) = u(AC)u f (B) 



(36) 



(37) 
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Proof: As the KMS state LOf is the bounded perturbation of the unique KMS state 

*(f) '(f) (f) 

lo and o~s commutes with oq for all t,s G R, ujf is a\ -invarinat and, thus, Lf is 

well-defined (cf. Corollary 2.3.17 of Ref. 0|). Then, because of the assumption on the 

spectrum of Lf, one has, for t — > oo, 

(^vr^aP^))^)-^,^/)^/,^^)^) = (^,e* L ^ B ) = J^JXe^^l^iA , 

where ip = ip — (Qf,ip)Qf and cJ5b = 7r /(-^)^/ — (Qf,TTf(B)0,f){lf are absolutely 
continuous vectors, {-E(A)} is the spectral family of Lf and the Riemann-Lebesgue 
lemma is used. Corollary 5.3.9 of Ref. |3| asserts that ffy is separating for the bi- 
commutant iTf{T)" of 7Tf(^ r ) since w/ is a KMS state. Then, as in the proof of 
Theorem 4.3.23 of Ref. 0], one concludes 

Km(iP,ir f (al f \B))(p) = (V, <p)(n f , ir f (B)Sl f ) = (il>,<p)u f (B) , 

from the above formula. As u is a bounded perturbation of u)f, it is expressed by a 
vector n G Hf. u{A) = (O , 7r/(^)fi )/||^o|| 2 (cf. Corollary 5.4.5 of Ref. 0) and, 
thus, one obtains the desired result: 

lim^iaf'fBlC) = Hm (n ,<Kf(A)TTf(al f) {B))ir f (C)n )/\\n f 

t—>oc f^oo 

= u f (B)(no,nf(A)TTf(C)n )/\\n f =u(AC)u f (B) . 



From this lemma, one immediately obtains 

L 

,lim u(r%f to a^_ to ((W f - Y ^N X ) 1, + kV)^^* 



Uf 



Wf-Y, ^x) ® ii + «f) o;(r2 f to f ? rS* 



A=l 
L 



((W/-5>AJVA)®li)+o;/(/riO . (38) 



A=l 



Moreover, |u;^(/cV")| < |/c|||V| 
In short, we have shown 



ttm Q T = w/((W> - Y W^a) <8 li) - w(K,(W> - J] Ma^VaX) <8> li) + 0(|«| 



T ^ + °° A=l A=l 



Remind that a; and ujf are K^-perturbed KMS states of the product states pi®u>cc and pf^ujcc-, 
respectively, where ujqc is the reservoir grand canonical state (a KMS state at with respect to 
t^&^L). Then, the stability of KMS states (Theorem 5.4.4 of Ref. @|) leads to ]im K -^Uf(A <g> 
li) = Tr(pfA) and lim K ^ w(-4 <g) li) = Tr(pjA) ( v ^4 G see also Ref. 0) and, hence, 

L L 

hm hm /?Q T = Tty,(w> - Y ^a) ~ ^ Pl (u t0 (Wf - Y 

= -Tr^/lnp/) +Trpi(u to (lnpf)u$ ) = -Tr(p f lnp f ) + Tr(p to (lnp f - lnp to )) + Tr(pilnpi 
= - - S{pf\p to ) . 



14 



Proof of Proposition 6: The formula (|3*Tj) is still valid: 



N 

-*E 

3=1 



A A 



(39) 



Let be the solution of T 



the solution of F, 
AW = 



"00 





1, 



d p"0') 



1 d r (i) 

A ' apt 



irp')f t (1) (W# ® li + kV) and V 



be 



((Wj-! + p(t)AHQ ® li + kV) with 



has f?V«(^)f?> 



(Wt — Wq)/N, then, by comparing these equations with that of Ft, one 
cP(A), u(rP(A)) = W (a2(5f (A))) = 

" E ^a) ® li + «V) - (Wj - E Vx^x) 



A 

Hence, we have 



a: 



li - 



(&-E 



-IUJ 



A 

f (D 



(^-Erf)(fE t0 )fS 



2],— V Tj-i+*o 



i K r T J - 1 ^ ) 1 (^-E^i 1) )(ff- ) ))r 

J A 

( r T,- 1+ ,oC 1+t o(4U((^ " E^A) ® ll + 



-to 



(r-i^-to^.,^ {<4?- to (W - E ma*a) ® ii + «v) }r^._ 1+to f 



(0)_* 
Tj-i-to 



where T t = 1 (r i ). On the other hand, the states u>,Uj and the evolutions a). 3 ' (j = 1, 2, • • • N) 
satisfy the assumptions necessary for the Lemma, we have 



(?) 



u.n r l "; I.'-, , «)'.' + .,(^)r*,_ i 



-Tj-i— 1 rj_i+to ij— to 1 



;(o)* 



r 

j-i+to — T 3 -_i— to 



r 



(0) 



-,(0)* 



lim w r 

Z>_i-*+oo V 

w ( r (0 i 



-fj-i-to t j-i+*o Tj_i+to -Tj-i-to 

(A)n 



ujj(A) =oJj{A) , 



(40) 



n(0)* 



lim w ( r> . . f P f , (A)) r% , , 

1 Tj_ 2 +to T,_ 2 +to V T i-i-*o V V Tj-2+to 



fj-i-to ^i-2+*o Tj-a+to -T,_2-to 



(41) 



Therefore, by taking the limits Tj — > +oo , Ij-i — > +oo in this order, we get 
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lim lim l-iuUh-^^Wf M) - Ts^^f, ,) r *f- 1 

J A A 

= <Oj({Wj - E fJ>\N\) ® l! + «y) " Wj _i(f ((Wj - E fi\N\) ® l! + K y)f f >) 



As before, one has 



max 



x(i)\/r"C7)\r"0')* 



"0> 



-Erf)(C j )f* Uj * < M «o ii(*i -E^i 1J )(^)ii 



)l} < l«l IMI , 



A 

< 2|«| t H^ll (W,- - E>a^a) ® li 



"0> 

to 



^(^(T^-E^A^f]®!!) 



(42) 
(43) 

(44) 



~(i) ~(i) 
where u\ is the solution of Uq 



dt u t 



iv^\Wj-i + (p(t)AW). More- 
over, the state LOj is a K^-perturbed KMS state of the product state pj <g> ojgc-, where 
Pj = exp(— (3(Wj — J2x MA^Y\))/"j with Ej the grand partition function, and, thus, lim K ^o Uj(A <g) 
li) = Th:{Apj) A € JTj). As a consequence, one gets 



lim lim lim 

k^O Tj - 1 ^oo Tj ^oo 



— IUJ 



(Si -Erf)(rf,)n ) +«, ( (*i - Erfj^jn 



:0> 
'to 



= Tr( Pj (^ - E^a)) - Tr^iuliHWj -^» x N x )i 

x x 
= /T 1 {-Tr( Pj \n Pj ) + Tr(p j _ 1 \np j )} = /T 1 {5( Pj ) - S^) - S{ P] \p 3 ^)} . 

where = pj^±u to , and, thus, 

AT 

lim lim lim ••• lim /3Qt = S(pf) - S(pi) S(pj\pj-i) , 

where P j = P o and pj = pN- 

Now we investigate the contribution from the relative entropies. Reminding that N\ commutes 
with Wo and Wf, one finds 

Stelfc-i) = /3<AugVug> >,-_! + /3((ugV J _ 1 ^> >,-_! - (W>_i>,-_i) + ln(eW- e -W) J _ 1 , 
where (-A)j-i = Tr(pj_iA) (\M G ^5). On the other hand, we have 

tigWug)* = e'^AWe-^- 1 * 1 + e iW ^" AW Aug** + An|f AW^f* , 



u 



— iWj-is 



W 3-Mo* ~ w j-i = i I ° ds(p{s)e lW i- lS [AW, Wj-^e 

t 

+i J° ds^(s) {e^- lS [AW, W^} A*#> + Au^ [AW, Wj-^uW*} , 



dTe TW ^AWe~ TW ^ 
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+ f dT [ T dr'e TW i- 1 AWe-<> T - T 'm- 1 AWe- T ' w f , 
Jo Jo 

where Au[^ = i / *° dsip(s)u^ AWe tW i- 1 ( t °~ s \ and, since pj-i commutes with Wj-i and ||Wj|| < 
\\Wo\\+j/N\\W f -Wo\\<\\W \\ + \\W f -W \\=K, 



7 Cj> 



< 2||AIU|| 2 



JO 



(vgWj.ivg* ) j - 1 -{W j - 1 ) j J <AK\\AW\\ 2 f°ds [ S ds'^sMs')] , 

JO " 

|( e W--i e -^) J -_ 1 -l + ^AW) J -_i| < ||AIU|| 2 I dr f dT'e Ki - T+T ' + \ T - T 'V , 



where we have used ([AW, Wj-i\)j-\ = 0. Therefore, if N is large enough, we obtain \S(pj\pj-\)\ < 
if'||AIU|| 2 = K'\\Wf - Wq\\ 2 /N 2 with some positive constant K', and the desired result: 



lim lim lim • • • lim (3Q T - {S(p f ) - S( Pi )} 



< 



K'\\Wf - W Q f 
N 



V. CONCLUSION 

In the first half of this article, we have shown that, if the evolution is L 1 -asymptotic abelian and 
the M0ller morphism relating a nonequilibrium steady state to a local equilibrium state is invertible, 
the natural non-equilibrium steady states can be regarded as MacLennan-Zubarev nonequilibrium 
ensembles. At first sight, invertibility of the M0ller morphisms seems to be too strong since one 
can easily find a counter example. However, for several systems, the division of the whole system 
without a finite part (i.e. the case where J-s = 0) provides invertible M0ller morphisms and, thus, 
we believe that Proposition 2 holds generically provided that the system is divided appropriately. 
Also it should be emphasized that Proposition 2 does not exclude the possibility that the class of 
natural nonequilibrium steady states is wider than that of MacLennan-Zubarev ensembles since 
the generator of the automorphisims defining the steady states takes MacLennan-Zubarev form 
only in a subset of its domain. 

In the second half, we have shown that the small system coupled with a single reservoir would 
follow 'thermodynamic' processes and satisfies the Clausius inequality/equality However, it should 
not be regarded as a dynamical proof of the second law of thermodynamics since we start from 
canonical ensembles which are very outcome of the second law. We think that the importance of 
this observation lies in the facts that (i) the dynamical evolution is consistent with thermodynamics, 
(ii) one can define thermodynamic heat microscopically, (iii) entropy generation is identified for 
step-wise evolution, and (iv) a characterization of a quasistatic change is given. 
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